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Meson transition amplitudes of the vector, axial-vector, and tensor quark currents are 
analyzed within dispersion formulation of the constituent quark model. The form factors in 
the decay region are given by relativistic double spectral representations through the wave 
functions of the inital and final mesons. We perform heavy quark expansion of the quark- 
model mesonic matrix elements with a next-to-leading order accuracy and demonstrate that 
matching this expansion to the heavy-quark expansion in QCD requires subtractions in the 
double spectral representations for the form factors and allows fixing the subtraction terms. 
The Isgur-Wise function and next-to-leading order universal form factors are calculated. 


I. INTRODUCTION 

Theoretical description of hadronic amplitudes of the quark currents is one of the key problems of particle 
physics as such amplitudes provide a bridge between QCD formulated in the language of quarks and gluons and 
observable phenomena which deal with hadrons. In particular, the knowledge of such amplitudes is necessary 
for extraction of the parameters of the quark-mixing matrix in the Standard Model from the experiments on 
weak heavy hadron decays. The difficulty in such calculations lie in the fact that hadron formation occurs 
at large distances where perturbative QCD methods are not applicable and nonperturbative consideration is 
necessary. 

Various theoretical frameworks more or less directly related to QCD have been applied to the description of 
meson transition form factors: among them constituent quark models jl] [llj, lattice QCD simulations Jr] |l8| , 
QCD sum rules |25| , analytical constraints j2(i|,^7j . 

Lattice QCD simulations is the most direct QCD based nonperturbative approach and thus should in principle 
provide most reliable results. However it faces serious technical problems with placing heavy particles on the 
lattice and performing calculations with inclusion of quark loops. So far direct calculations with 6-quark are 
not possible and extrapolation in the quark mass is used which considerably reduces the predictive power of 
the method. 

Various versions of QCD sum rules to meson decays give very uncertain predictions strongly dependent on 
the technical subtleties of the particular version. A recent analysis disregards the three-point QCD sum 
rules in favor of the light-cone sun rules which however involve more phenomenological inputs. 

Constituent quark models have proved to be a fruitful phenomenological method for the description of heavy 
meson transition form factors. Constituent quark models employ the fact that heavy meson consists of a 
heavy quark and light degrees of freedom with have quantum numbers of the quark state and assume that 
these light states can be approximately described by an effective constituent quark. The application of the 
various versions of the constituent quark picture to the decay processes has a long history. The first models 
were based on a semirelativistic |lj or nonrelativistic [^] Q considerations and didn’t treat the quark spins 
thoroughly. Further developments demonstrated the relativistic description of the quark spins to be important 
for consistent applications to meson decays. The exact solution to a complicated dynamical problem of treating 
the spins of the interacting particles is not known, but a simplified self-consistent relativistic consideration 
of the quark spins can be performed within the light-cone quark model (LCQM) [(i||. The difficulty with the 
application of this approach to the decay processes lies in the fact that the applicability of the model is restricted 
by the condition q 2 < 0, while the physical region for hadron decays is 0 < q 2 < (Mi — Mf) 2 , Mij being the 
initial and final hadron masses, respectively. The problem is connected with the contribution of the so-called 
pair-creation subprocesses which cannot be taken into account thoroughly in the model, except for a trivial 
case of a pointlike interaction. At spacelike momentum transfers the contribution of these subprocesses can 
be killed by choosing an appropriate reference frame, whereas at timelike momentum transfers such a frame 
does not exist and thus pair creation contributes together with the partonic contribution. An estimate of the 
form factors at timelike momentum transfers within the LCQM has been done in Q where the transition form 
factors have been calculated in the spacelike region and then numerically extrapolated to the timelike region 
assuming some particular q 2 behavior. This works well if the accessible momentum transfer interval is not 
large. However, in the heavy-to-light meson decays this interval is of order of the heavy meson mass squared 
and the extrapolation procedure yields large uncertainties. The approach of refs [|ll] [l3| calculates the partonic 
contribution to the form factors at timelike momentum transfers and thus avoids extrapolations. However, the 
nonpartonic contribution is omitted. Such treatment can be justified if the latter is small that is not obvious 
a priori. Thus, one can see that for a reliable description of the decay processes within the LCQM it might be 
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reasonable to find another formulation of the model appropriate also at timelike momentum transfers. Such a 
formulation has been proposed in 0- 

The approach of |l4|] is based on the dispersion formulation of the LCQM. Namely, the transition form factors 
obtained within the light-cone formalism at q 2 < 0 iQ, are represented as dispersion integrals over initial and 
final hadron masses through their light-cone wave functions. The spectral densities of these representations can 
be calculated from the Feynman graphs. The transition form factors at q 2 > 0 are derived by performing the 
analytical continuation in q 2 from the region q 2 < 0. As a result, for a decay caused by the weak transition of 
the quark Q{rrii) — > Q(m/), form factors in the region q 2 < (m» — w/) 2 are expressed through the light-cone 
wave functions of the initial and final hadrons and can be calculated in the decay region avoiding dangerous 
extrapolations. 

At the same time, such a dispersion formulation of the LCQM allows also a bit different view on the problem 
of constructing the form factors within the quark model: actually, we calculate the double spectral densities 
from the Feynman graphs and thus obtain unsubtracted spectral representations. On the other hand we know 
that in general subtractions might be added to such representations and we need some additional arguments to 
decide on the necessity and the structure of the subtraction terms. 

Let us notice that once we are working within an approach not directly deduced from QCD it is important 
to preserve essential features of the underlying fundamental theory in the model. Thus matching the results 
obtained within the quark model to rigorous QCD results might be helpful for bringing more realistic features to 
the model. We use such matching to QCD results for determining the subtraction terms in the double spectral 
representations of the form factors. 

There are few QCD predictions on meson form factors in the decay region. Fortunately, there is an important 
case in which QCD provides rigorous and model-independent results on transition form factors: namely, meson 
decays induced by transitions between heavy quarks. In this case the heavy-quark symmetry simplifies all 
the form factors in the leading 1/toq order to the universal Isgur-Wise function (^8|. A systematic expansion 
in inverse powers of the heavy quark mass can be constructed within the Heavy quark effective theory [0 
based on QCD with the heavy quarks. Although the form factors appearing in different 1/toq orders can not 
be determined within the HQET, the latter restricts the structure of the expansion and the number of these 
universal form factors in each 1 /ttiq order and provides constraints on some of them at zero recoil. 

The idea of matching the quark model form factors to the structure of the heavy-quark expansion of HQET has 
been employed in the ISGW2 model ||. However the ISGW2 model does not calculate the form factors within a 
relativistic dynamical approach on the basis of considering the meson structure but rather gives prescription for 
constructing the form factors. The approach of ref. |Q also employs matching to the heavy-quark expansion and 
obtains interesting results within the constutuent quark picture assuming a strong peaking of the momentum 
distribution of the constituent quarks inside a meson and avoiding details of the dynamics of the process. This 
model imposes constraints on the form factors but does not allow one to calculate them. 

The idea implemented in this work is to perform the heavy-quark expansion in the dispersion formulation 
of the quark model and to match this expansion to the HQET. This matching is used for determining the 
structure of the subtraction terms in the spectral representations. We consider the form factors of pseudoscalar 
meson transitions into pseudoscalar and vector final mesons through the vector, axial-vector, and currents and 
analyze the expansion in the leading and next-to-leading order. Further constraints on the subtraction terms 
are obtained by considering the case of the heavy-to-light quark transition and analyzing the relations between 
the form factors of the tensor, vector and axial-vector currents. Comparing the relations obtained within the 
dispersion quark model to the general relations found by Isgur and Wise 1321 allows us to fix the subtraction 
terms. 

Our main results are the following: 

1. Applying the approach of |29| to heavy meson transition throught the tensor current, we derive the de¬ 
composition of the form factors h g+ , h g _, h go , and h s (the definitions are given in the next section) in the 
next-to-leading order in 1/toq. In particular, we have found that the 1/toq corrections to the form factor h g+ 
at zero recoil vanish just as for the form factors hf + and hf (the Luke theorem). Although the heavy meson 
transitions through the tensor current do not correspond to a real experimental situation we need the expansion 
for such form factors as a bench mark for comparison of the quark model. 

2. We present an improved dispersion formulation of the LCQM: the transition form factors of pseudoscalar (P) 
mesons into pseudoscalar (P) and vector ( V ) mesons induced by the vector, axial-vector and tensor currents 
are given by double spectral representations through the wave functions of the initial and final mesons. The 
unsubtracted double spectral densities are calculated from the Feynman graphs and the subtraction terms in 
the spectral representations are determined by matching the 1 /toq expansion in the quark model to the HQET. 

3. We perform the heavy-quark expansion of the form factors of the dispersion quark model up to NLO and 
calculate the Isgur-Wise function and the NLO universal form factors assuming a strong peaking of the soft 
wave function near the qq threshold with a width of order A qcd- We obtain in the LCQM the relations 


y 2 (w) = 0; X3(w) = 0; 6(w) > 0, £ 3 (1) = (z)/ 3, 


( 1 ) 
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where (z) is an average kinetic energy of the light quark in the heavy meson rest frame. 

We observe that matching the 1 /toq expansion of the quark model to HQET in LO and NLO requires 
subtractions in some of the form factors describing the P —> V transition and restricts the form of the subtraction 
terms. 

4. We consider the heavy-to-light meson transitions in which case a small parameter Aqcd/wq emerges 
and analyze the form factors in the leading A qcd/i^q order. Requiring the fulfillment of the Isgur-Wise 
relatins for the heavy-to-light transitions [ |32| further constrains the subtraction terms providing explicit spectral 
representations with subtractions for the form factors of interest. 

In the next section we remind the definitions of the form factors and list the HQET results for the amplitudes 
of the P —> P and P —> V transitions through the vector, axial-vector and tensor currents up to NLO. In 
section 3 we present the quark-model results for the meson decay form factors at arbitrary masses: namely, we 
calculate the unsubtracted spectral densities and explicit form of the subtraction terms which form is motivated 
in the subsequent section. In Section 4 the 1 /toq expansions of the wave functions and the form factors in 
quark model for the heavy-to-light transitions are considered. In Section 5 the heavy-to-light case is discussed. 
Section 6 illustrates the main results with numerical estimates and evaluate the Isgur-Wise function for various 
quark model parameters. Conclusion summarizes the results. 


II. MESON TRANSITION AMPLITUDES AND HEAVY-QUARK EXPANSION IN QCD 


The amplitudes of meson decays induced by the quark transition g 2 —> qi through the vector V g = 917 g q 2 , 
axial-vector A g = < 717 ^ 75 ( 72 , and tensor = qicr gv q 2 currents have the following structure 1 32 1 


< P(M 2 ,p 2 )|^(0)|P(Mi,pi) > = fi(q 2 )pi g + / 2 (g 2 )P 2 M , 

< V(M 2 ,P 2 , e)\Vf l (0)\P(Mi,pi) > = 2 g(q 2 )e fJ , vaf) e*''p^ p%, 

< V(M 2 ,P2,e)\A^(0)\P(M 1 ,p 1 ) > = ie* a [ f(q 2 )g,j, a + ai(q 2 )p la px^ + a 2 (q 2 )pi a P2^}, 

< P(M 2 ,p 2 )|T„„( 0 )|P(Mi,pi) > = - 2 is(q 2 ) {pi g p 2v ~Pi„P 2 »), 

< V(M 2 ,P2,e)\T tlv (0)\P(M 1 ,p 1 ) > = ie* a [gi{q 2 )e gva 0P 10 + g 2 (,q 2 )e^app 213 + 9o(q 2 )pi a e^p-yPiP2 ], (2) 

with q = pi — p 2 , P = pi + p 2 - We use the notations: 7 s = * 7 ° 7 1 7 2 7 3 , a gu = §[ 7 ^, 71 ,], e 0123 = — 1 , and 
Sp( 7 5 7 M7^7« 7 /3) = 4 ie^ a P. 

The relativistic-invariant form factors contain the dynamical information on the process and should be cal¬ 
culated within a nonperturbative approach for any particular initial and final mesons. 

For analysing the transition in the case when both the parent and the daughter quarks inducing the meson 
transition are heavy, i.e. mi — m 2 A qcd it is convenient to introduce a new dimensionless variable 

to = V 1 V 2 = Ul j u^m 2 (/ anc ^ velocity-dependent form factors connected with 4-velocities and not 4-momenta as 
in (|j) in the following way 


< P(M 2 ,P2)|^(0)|P(Mi,pi) > 
< E(M 2 , P 2,e)|H M (0)|P(Mi,pi) > 
<y(M 2 ,p 2 ,e)|A^(0)|P(M 1 ,p 1 ) > 
< P(M 2 ,P2)\T^(0)\P(M 1 ,p 1 ) > 
< V(M2,p2,e)\T^(0)\P(M 1 ,pi) > 


= \]MiM 2 [hf + (u>)(v 1 + v 2 )v + hf_ (uj)(vi - 
= VMiM 2 h g (ui)e tlvct ge* v v ? v$, 

= ie* a \/MiM 2 [hf(u)(l + u)g lia - h ai {uj)vi a v lfJ , - h a2 (u!)vi a v 2g .], 
= - 2 i \/M 1 M 2 h s (uJ ) {vi g V2v - Vl v V2n), 

= it a \J MiM 2 [hg + (u)e^ a f 3 (vi + v 2 Y + h g _ (uj)e gva0 (vi - v 2 ) /3 
+h go (w)v lae^g-yvfv^]. 


These form factors are related to the form factors introduced by the relations (E) as follows 


(3) 


fi = 

h = 




VM 1 M 2 

Mi 


[hf + +h f _], gi = 


M 2 


\/ M 1 M 2 

g = - r l 1 -= h g , s = 


[hf + -h f _], g 2 = 

1 


2 y/ AI\ M 2 


2y/ M 1 M 2 


\J M\ M 2 

Mi 
ho , 


[hg + + hg _\ 7 

[hg+ ~ h g _], 


f = \J M]M 2 (1 +u>)hf, a 2 = — 


1 


\/ .V/ [ A / 2 


ha .2 ) ^1 — 


1 


M 2 


1 


1 


y/MiMi Ml 


9 0 = 


yJM x M 2 Ml 


i 9o 


(4) 
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The form factors h are convenient quantities as in the leading 1 /mg order all of them are expressed through a 
single universal function of the dimensionless variable us - the Isgur-Wise function. A consistent heavy-quark 
expansion of the form factors, i.e. expansion in inverse powers of the heavy-quark mass, can be constructed 
within the Heavy Quark Effective Theory based on QCD with heavy quarks. 

The general structure of the 1/toq expansion of the heavy quark form factors in QCD for the meson 
transition M\ —> M 2 induced by heavy quark transition m 2 —■> m\ have the form (omitting corrections 
0(a s , a s /m Q , 1 /m 2 Q ): 

h f+= £ + 
hf _ = 

hg = £ + 

hf= £ + 

ha 2 = £ + 
h s = £ + 
h g+ = — £ - 

V = 


In the leading 1 /toq order (LO) all the form factors are represented through the single universal Isgur-Wise 
function £(w) , whereas in the next-to-leading order (NLO) the 4 new form factors pi,p 2 ,X 2 , and £3 appear. 
The universal form factors are functions of a single variable us. 

The form factor £3 originates from the expansion of the transition quark current, and the form factors pi, p 2 , X 2 
are connected with the nontrivial relationship between the mesonic states in the full and the effective theory. 
The universal form factors satisfy the conditions 

£(!) = !, Pi(l)=P2(l) = 0, (6) 

whereas no constraints on £3 and \2 are imposed by the heavy quark symmetry. As found by Luke p9| , the 
1 /toq corrections to the form factors hf and hf + vanish due to kinematical or dynamical reasons. One can see 
that the same is true for the form factor h g+ : namely, 

h g+ (l) = -l + 0(l/m 2 Q ). (7) 

The parameter A in (j^) comes from the 1 /toq expansion of the mass of a meson consisting of the heavy quark 
and light degrees of freedom 

Mq = m Q + A—I- 0(l/m Q ). (8) 

In our notations for heavy quarks and mesons, this gives 

Mi = m 2 + A + ..., M 2 = mi + A + ..., (9) 

for the parent and daughter particles, respectively. 

It is straightforward to derive the following useful relations 

Mi + M 2 _ m 2 + mi 
sJMiM 2 A /miW 2 


" / 1 

i \ 

( m 2 - TOi N 

\ A 


1—1 — + 

— 

— 

) 7T • 


\m i 

m 2 / 

\ TO2 + TOi j 

/ 2 



1 1 

mi m 2 

1 1 

mi m 2 

1 1 

TOi 7712 

1 1 


Pi, 




A^ i 1 . 

77? 4- P 2 H- (P 1 - £ 3 ) 

2 mi 


mi m 2 

1 1 


m 2 


u-lA 1 1 

——T7t£ d- P2 H- 

us +1 2 mi m 2 


Pi 


us — 1 
us + 1 


£3 


TO 1 US 

us — 1 1 


Y (-M + 2 (us + i)x 2 - 6), 

1 


US + 1 TOi 

1 1 


TOi 


m 2 

A 


A J_ 

2 TOi 


P2 - 2X2 - 


1 


US + 1 


& ) H- (pi - 6)- 


m 2 


1 


m 2 


-Pi 


m 2 

1 

mi' 


r£ - C 3 + Pi 


-P2, 


_l _ l_\ A J_ 

TOi TO 2 J 2 + 777-2 3 ’ 


TOi \ US + 1 


2X2 
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Mi - M 2 
\J MiM 2 

\J Mi M2 
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i/mim 2 


y'rnS 
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1 > 
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' 

1-1 — 

H- 
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\m 1 
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1 > 

A 

1 + ( — 

H- 
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(10) 


where the dots denote higher order terms. 

Using the relations (|]) and (jl(]), we obtain for the form factors i) the following expansions 
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A 

s 


a 1 


92 
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CLi 


a 2 


/ 

go 


mi 


V m i m 2 
m 2 


y/mim 2 

1 

2 i/rrTpm^ 
mi 


£ H-(pi + £ 3 ) H-(pi - £ 3 ) 

mi m 2 

£ H-(pi - £ 3 ) H- {pi + £ 3 ) 

mi m 2 

£ H-(Pi - &) 4-(pi - 6) 

mi m 2 


VmS 

m 2 

yjmim 2 
1 

2 yjmim 2 
1 


£■ 


1 1 / 

—P 2 H-(Pi ~ £ 3 ) 

mi m 2 


£ H-P 2 H-(pi + £ 3 ) 

m 1 m 2 

£ H-P 2 H-(pi - £ 3 ) 

mi m 2 


1 1 


v /mim 2 m 2 w + 1 

1 


[—A£ + 2 (w + 1)X2 - £3] 


i/mim 2 

y/mim 2 (uj + 1 ) 


1 , . , 1 , _ A 

— (Pi - £3) H-(Pi - £3)- — 

mi m 2 mi u> + 1 


2x2 + 1 u 

mi mi u> + 1 


£3 


£ + 


A/ 1 


u> + 1 \ m i 


1 

m 2 


1 1 
— P 2 H- 


mi 


m 2 


Pi 


uj — 1 

w + 1 


£3 


1 


(mim 2 ) 3 / 2 u + 1 


A£ + £3 


2x2 


( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


For the following analysis it is worth noting that the behavior of the combination 2pip 2 ■ g — mi ■ g 2 and / in 
LO and NLO coincide, namely 


/ ~ 2pip 2 ■ g - mi ■ g 2 . 

It is also convenient to introduce the form factor a' 2 such that 

CL 2 = cl 2 2s. 

In what follows we need the expansions of the following linear combinations of the form factors a' 2 and ai 

, 1 


aim 2 — a 2 m 1 = — 
ai m 2 + a 2 mi = — 


y/mim 2 

1 

yjmim 2 


[4X2 - £ 3 ] 


2A£ w - 1 
£3 


U) + 1 


UJ + 1 


( 21 ) 

( 22 ) 

(23) 

(24) 


III. TRANSITION FORM FACTORS IN THE DISPERSION QUARK MODEL 

The results presented in the previous section are strict consequences of QCD in the heavy-quark limit which 
however cannot provide more information on the universal form factors £ and Pi,P 2 ,X 2 )£ 3 - They must be 
calculated within a nonperturbative dynamical approach. We study the form factors within the dispersion 
formulation of the quark model [ ffdlfTAI which has proved to be a reasonable framework for describing meson 
decays. We start with q 2 < 0 and represent the form factors as double spectral representations in the invariant 
masses of the initial and final qq pairs. The form factors at q 2 > 0 are derived by performing the analytical 
continuation. 


5 





























































The transition of the initial meson q{rn 2 )q{m ^) with the mass M\ to the final meson q(rni)q{mz) with the 
mass M 2 induced by the quark transition m 2 —► m\ through the current g(mi) J M g(m 2 ) is described by the 
diagram of Fig.jl]. For constructing the double spectral representation we must consider a double-cut graph 
where all intermediate particles go on mass shell but the initial and final mesons have the off-shell momenta p\ 
and and p 2 such that p\ = Si and p 2 = S 2 with (pi — p 2 ) 2 = S 3 = q 2 kept fixed. 


/ 



FIG. 1. One-loop graph for a meson decay. 


For the transition B —> D,D* m 2 = m&,mi = m c , and m 3 = m u . The constituent quark structure of the 
initial and final mesons is given in terms of the vertices Fi and r 2 , respectively. The initial pseudoscalar meson 
vertex has the spinorial structure Ti = ijsGi/y/Nc; the final meson vertex has the structure T 2 = /ysG^/v^Ve 
for a pseudoscalar state and the structure r 2/i = [A r y ll + B{k\ — ^ 3 )^,] G^/y^/Vc, A = — 1 , B = l/(y / S 2 +mi + 7713 ) 
for an S-~wa,ve vector meson. 

The double spectral densities / of the form factors are obtained by calculating the relevant traces and isolating 
the Lorentz structures depending on p\ and p 2 . The invariant factors of such Lorentz structures provide the 
double spectral densities / corresponding to taking into account contributions of the two-particle singularities 
in the Feynman graph. Let us point out that this procedure allows one to obtain unsubtracted double spectral 
densities, whereas subtraction terms should be determined independently. In this paper we determine the 
subtraction terms from matching the 1 /mg expanded form factors of the quark model to the heavy-quark 
expansion in QCD. 

At q 2 < 0 the spectral representations of the form factors have the form 0 


+ / 2 \ 

oo sj(s 2 ,q ) 

(mi+m 3 ) 2 s~(s 2 ,q 2 ) 


(25) 


where the wave function = Gi(si)/(si — M 2 ) and 


si{s 2l q 2 ) = 


s 2 (m 2 +m 2 — q 2 ) + q 2 {m\ + m 2 ) — (mf — m 2 )(m 2 — m 2 ) A 1 / 2 (s 2; m 3 , m\)\ 1 ^ 2 (q 2 , m \, m 2 ) 


2m 2 


2 mf 


and A(si, s 2 , S3) = (si + S2 — S3) 2 — 4 sis 2 is the triangle function. 

The unsubtracted double spectral densities //(si, s 2 , g 2 ) of the form factors read [ p~i||lE| ]: 

s = 2 [mi 0:2 + m. 2 a 1 + m3(l — aq — 0 : 2 )], 
fi = 2mis + 4a 2 [s2 - (mi - m 3 ) 2 ] - 2m 3 s, 
h = 2 m 2 s + 4«i[si - (m 2 - m 3 ) 2 ] - 2m 3 s, 
g = As — ^Bfj, 

gi = Afi - 8/3 + 82? (mi + m 3 )/?, 
h = Af 2 + 8 B{m 2 - 7713 )/?, 
d 2 D = — 2s + 4BC 2 a 1 + ctq 2 Co, 
am = —4 A (2m 3 + BC{)at\ + an Cb, 

f D = - 4 A[mim 2 m 3 + -y-(s 2 - m 2 - m 2 ) + ^y-(si - m 2 - m 2 ) - ^(s 3 - m 2 - m\)\ + Co/ 3 , 
goD = -8A«12 - 8 B [—m 3 aq + (m 3 - m 2 )an + (m 3 + mi)ai 2 ], 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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where 


ai = [(si + s 2 - s 3 )(s 2 - ml + m 2 ) - 2s 2 (si - m\ + m§)] /A(si, s 2 , s 3 ), (36) 

o<2 = [(si + s 2 - s 3 )(si - m\ + mg) - 2si(s 2 - m\ + mg)] /A(si, s 2 , s 3 ), (37) 

f3 = ^ [2mg - ai(si - mj + mg) - a 2 (s 2 - to 2 + mg)] , (38) 

an = a? + 4/3 s 2 /A(si, s 2 , s 3 ), a i2 = a 3 a 2 - 2/3(si + s 2 - s 3 )/A(si, s 2 , s 3 ), (39) 

Co =-8A(m 2 - m 3 ) + 4BC 3 , Ci = s 2 - (m 3 + m 3 ) 2 , (40) 

C 2 = si - (m 2 - m 3 ) 2 , C 3 = s 3 - (mi + m 2 ) 2 - Ci - C 2 . 


We label with a subscript ’D’ the double spectral densities of the form factors which require subtractions. 
We fix this subtraction procedure by requiring the 1 /ttiq expansion of the form factors to have a proper form 
in accordance with QCD in leading and next-to-leading orders for the case of a meson transition caused by 
heavy-to-heavy quark transition. As we shall see this yields the double spectral densities which include properly 
defined subtraction terms 


/ — Id + [(Mf — Si) + (M 2 — s 2 )]<7, 


hi = diD + 


d 2 = a-in + 


go = goD + 


(u) + 1 )m 2 V v 7 ®! 


Mi — si | Mf — s 2 \ g 
2 ’ 


V®2 


1 /Mi — si M| -s 2 \ j 


(w + l)mi V v 7 ®! 


V 7 ®^ 


2’ 


(w + l)mim 2 \ V ^T 


Mf - si ( M| - s 2 A g_ 
v 7 ^ 7 2' 


(41) 

(42) 

(43) 

(44) 


As the analytical continuation to the timelike region is performed, in addition to the normal contribution 
which is just the expression (|2^) taken at q 2 > 0, anomalous contribution emerges ||T9[. The origin of the 
anomalous contribution is connected with the motion of a zero of the triangle function A(si,s 2 ,s 3 ) from the 
unphysical sheet at q 2 < 0 onto the physical sheet at q 2 >0 through the normal cut between the points 
sr(®2,® 3 ) and S+(s 2 ,S 3 ). 

To be more specific, we can write 


A(si, s 2 ,s 3 ) (si Si)(si Si ), 


(45) 


where sf = (a/s^ — y 7 ®^) 2 an< ^ ®f = (v 7 ®^ + y 7 ®^) 2 - At q 2 < 0 both singularities sf and sf are located on 
the unphysical sheet and do not contribute to the spectral representation of the form factor. As q 2 becomes 
positive, sf moves around s(" and appears on the physical sheet. The pinching of the points sf and sj~ occurs 
at s 2 such that sf (s 2 ) = sj~(s 2 ). The corresponding value s 2 (s 3 ) reads 


s 3 + m 2 — mf 
2 a/s 3 


i 2 9 A ^ 

S 3 + mf — 777-2 


2^/53 


+ (to§-to 2 ), s 3 <(m 2 -mi) 2 . 


(46) 


The spectral representation of the form factor at q 2 > 0 takes the form U4 


OO °1 

f(Q 2 ) = J V2{S2) J <Msi) 


(mi+m 3 ) 2 


f{s 1 ,s 2 ,q 2 ) 
A 1 / 2 (si, s 2 ,g 2 ) 


OO *1 _ 

hII (47) 


One should also take into account that if the double spectral density /,; is singular at the point sf then a 
properly defined spectral representation contains also another kind of subtraction terms: these terms appear as 
the contribution of a small circle around sf which is the lower integration limit in the anomalous part. These 
subtractions have quite different nature and different form than subtractions in the spectral representations at 
q 2 < 0. The corresponding expression can be found in ]l5| . 

The normal contribution dominates the form factor at small timelike and vanishes as q 2 = (m 2 — mi) 2 while 
the anomalous contribution is negligible at small q 2 and steeply rises as q 2 —*■ (m 2 — mi) 2 . 

For pseudoscalar and vector mesons with the mass M built up of the constituent quarks m q and m qi the 
function p is normalized as follows ||l^| 


1 

87T 2 




(■ m.q - m q ) 2 } = 1. 


(48) 
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This equation is the normalization of the elastic charge form factor at q 2 = 0. 
The meson wave function can be written in the form 


V?(s) 


j,^es=3E ' W ( t > 

\[2 yjs (mi - in 2 ) 2 s 3 / 


k = 


A 1 /^ 


s, mj, m : 


2i/s 


(49) 


where w(k ) is the ground-state 5-wave radial wave function. 

In the next sections we analyze the form factors given by the dispersion representation ( ff?| ) with the spectral 
densities (^ 6 | pT| ) and (|4l|-[f4|) and demonstrate them to have the following properties in the case of a heavy 
parent meson: for the transition induced by the heavy-to-heavy quark transition they satisfy the LO and NLO 
relations p9| of the 1/niQ expansion in accordance with QCD provided the functions ipi are localized near the qq 
threshold with the width of order A qcd- For the meson decay induced by the heavy-to-light quark transition 
they satisfy the LO relations between the form factors of the vector, axial-vector and tensor currents |3^|. 


IV. HEAVY-QUARK EXPANSION IN QUARK MODEL FOR HEAVY-TO-HEAVY TRANSITIONS 

In this section we consider the form factors of the dispersion quark model in the case when both m 2 and 
nii are large. We calculate the universal form factors and demonstrate that requiring the structure of the 
1/niQ expansion in the quark model to be consistent with the structure of such expansion in QCD allows us to 
determine the subtraction terms. 


A. Soft wave function and normalization condition 

First, we need to specify the properties of the soft wave function of a heavy meson. A basic property of 
such soft wave function ip(s 1 mQ,mq, A) is a strong peaking near the qq threshold. For elaborating the 1/mq 
expansion, it is convenient to formulate such peaking in terms of the variable z such that s = (niQ + m 3 + z) 2 
(hereafter we denote the mass of the light quark as m 3 ). The region above the qq threshold which contributes 
to the spectral representation corresponds to z > 0. A localization of the soft wave function in terms of z means 
that the wave function is nonzero as z < Aqcd- In the heavy meson case itiq 3> m 3 ~ 2 ~ A. Let us notice 
that for a heavy meson the localization in terms of z is equivalent to the localization in terms of the relative 
momentum in the meson rest frame 


k 2 = z(z + 2ms) + 0(l/mQ). (50) 

The normalization condition ( p 8 | ) which is a consequence of the vector current conservation in the full theory 
provides an (infinite) chain of relations in the effective theory. Namely, expanding the soft wave function in 
1 /mq as follows 


ip(s,mQ,m 3 , A) = 


y/rnQ 


4>o{z, m 3 , A) 


m 3 


1 + ^m ^ Xl ^ Z,m3 ' A ) + 


we come to the normalization condition in the form 
[ dz(j>l(z)\fz(z + 2m 3 ) 3 ^ 2 


, , m 3 , , m 3 

1 + 7, -Xl(~) — -b • • • 

2 raq 2uiq 


= 1 . 


(51) 


(52) 


This exact relation is equivalent to an infinite chain of equations in different 1 /mq orders. Lowest order relations 
take the form 


/ dz 4>l(z)y/z(z + 2 m 3 ) 3//2 = 1, (53) 

/dz </>o {z)y/z(z + 2m 3 ) 3/2 xi(z) = 1. (54) 


B. The variables ui and a 1 

In the description of the transition processes the dispersion formulation of the quark model has the following 
feature: since the underlying process is the quark transition, the relevant kinematical variable for the description 
of the dynamics of the process is the quark recoil Co which is defined as follows 














q 2 = (m2 — mi) 2 — 2mim 2 (w — 1). 


( 55 ) 


The relationship between to and to is given by the condition that the spectral representation for the form factor 
is written at fixed value of q 2 and hence 

q 2 = (Mi — M 2 ) 2 — 2MiM 2 (to — 1) = (m 2 — mi) 2 — 2mim 2 (w — 1). (56) 

In the case of heavy particle transitions these quantities are related to each other as 

to = to + A (- 1 -J [to — 1) + O(l/mo). (57) 

\mi m 2 / 


We shall obtain the representations of the form factors as functions of the variable Co. The variables io and 
Co are different by the terms of order 1 /mg at nonzero recoil. On the other hand, the quark and meson zero 
recoil points coincide with 1 /mg accuracy. This means that in the analyses of the 1 /mg expansion at nonzero 
recoil the difference between the to and to might be ignored in the Isgur-Wise function, but gives nontrivial 
contribution to the NLO form factors. At the same time, at zero recoil the difference between to and to might 
be neglected both in the leading and next-to-leading orders. Namely, the quark model provides the expansion 
of the form factor in the following form 


h = ho(to) H— —h^fto) H— —hY’(to) ■ 

mi “~- 


1 fo 2 ) I 
m 2 


ho(to) + h' 0 (oo )A (— + — ^ (to — 1 ) + — h^\to) + —h^\to) + — 
\mi m 2 / mi m 2 


(58) 

(59) 


As we shall see later, among the NLO form factors only pi 2 are affected by the the difference between to and 
Co whereas £,£ 3 , \2 are not. 


C. Relative magnitudes of the normal and the anomalous contributions 


We are going now to demonstrate that the anomalous contribution comes into the game only in close vicinity 
of the zero recoil point whereas beyond this region is negligible. 

Let us study the behavior of the anomalous contribution in the region 

Co — 1 ~ mg^ 2+s \ (60) 

Introducing the variables z\ and z 2 such that si = (?n 2 + m 3 + z 3 ) 2 and s 2 = (mi + m 3 + z 2 ) 2 we find that the 
magnitude of the anomalous contribution is controlled by the value of z 2 (Co) such that = (mi + m 3 + z 2 (Co)) 2 
which is the lower boundary of the Z 2 integration. If z 2 (Co) becomes large, i.e. of the order ijiq, the anomalouis 
contribution is suppressed by the fall-down of the soft wave function. This suppression is at least stronger than 
1 /m,Q. This means that the anomalous contribution is nonvanishing only if 


z° 2 (u>) 


mim 2 Vw 2 — 1 + mim 2 w — m 2 
\/ m 2 + m\ — 2 m. 1 m. 2 C 0 


mi + 0(m 3 ) ~ A. 


(61) 


In the region to — 1 ~ mQ 2+e \ one finds z 2 (to) ~ mq^ . Hence the anomalous contribution comes actually 
into the game only in the 0(l/mg) vicinity of the zero recoil point but otherwise might be neglected. On the 
other hand, at the quark zero recoil point Co = 1 , the normal contribution vanishes and the form factor is given 
by the anomalous contribution. 

We shall calculate the form factors in the region Co — 1 = 0(1) where only the normal contribution should be 
taken into account in leading and subleading orders. 


D. The LO analysis 


To perform the LO analysis of the form factors let us start with the integration measure. With 1/mQ accuracy 
it can be represented in the form 


1 dsids20 (s 2 > (mi + m 3 ) 2 ) 9(s 1 < si < sf) 
167 t 2 A 1 / 2 (s 1 , s 2 , s 3 ) 


~ -^dz 2 \/Z 2 (z 2 + 2m 3 )^-9(z 2 )0(p 2 < 1 ), 
47 T z Z 


(62) 
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and the expression for z\ reads 


Z\ = z 2 uj + m 3 (uj - 1) + r\\J z 2 (z 2 + 2m 3 )\/ uj 2 - 1 + 0(1/mo). 


(63) 


Let us point out that the LO integration measure is symmetric in z\ and z 2 . 

Next, we need expanding the spectral densities (Eq-Bq). To this end we must take into account that under 
the integral sign z\ and z 2 are localized in the region z < A due to the soft wave functions <j>(z). 

In LO the kinematical coefficients (|36| ^) in the region uj — 1 = 0(1) simplify to 

(64) 

(65) 


A(si, s 2 , q 2 ) = Am\m\(uj 2 - 1), 
1 


a i = 


a- 2 = 


m 2 (uj + 1) 

1 


mi(u + 1 ) 


Z 2 ~ Zi 

m 3 + z 2 + —-— 

uj — 1 

z\ — z 2 

m 3 + Zi + — -— 

uj — 1 




OL\\ — + 


UJ 


^—(m 3 + zi)(m 3 + z 2 ) + ("t ^ 

+ 1 uj 2, — 1 


P 


^ 2 (uj 2 ~iy 


ol i 2 = o;ia 2 — 


(3 uj 


C 0 = —4m 2 (uj — 1), B = 


2m i ’ 


m.\m 2 (uj 2 — 1) ’ 

Oi = 2m\Z 2 , C 2 = 2m 2 (zi + 2 m 3 ). 


( 66 ) 

(67) 

( 68 ) 

(69) 


One finds the LO behavior of the form factor densities 
The latter reads in the LO 


5 ~ 2 ( m 3 + 

The LO expression for / takes the form 


|) is determined by the term proportional to s. 
z\ + z 2 + 2m 3 ' 


f D = (uj + 1 ) ( m 3 + 


uj + 1 


Zi + z 2 + 2 m 3 
uj + 1 


(70) 


(71) 


The spectral densities di and gon vanish in the leading order. 

Hence the LO relations (jlT] are fulfilled with the Isgur-Wise (IW) function 

l 

= J dz 2 cj) 0 (z 2 )y/z 2 (z 2 + 2 to 3 ) J y^o^i) (m 3 + ' 


(72) 


In ( [72| ) we used the equality of uj and uj with 1 /rag accuracy. The normalization condition of the LO wave 
functions (|5^) yields £(1) = 1. For the slope of the IW function at zero recoil, p 2 = —£'(1), one finds 

p 2 = ^ J dzy/z(z +2m 3 ) 3/2 (4>o(z)) 2 z(z + 2m 3 ) (73) 

Let us point out that the subtraction terms in the spectral densities do not contribute in the LO relations. 
As we shall see later, they are important in the NLO analysis. 


E. The NLO analysis of the form factors s, /1, f 2 , g 1 , g 2 , and g 


First, let us concentrate on the NLO relations (jlT] |l 6 ]). It is convenient to analyze the linear combinations of 
the form factors which do not contain the LO contribution. These combinations are 


g-s 

91 + fi 

92 + f 2 
fi ~ 2 m-is 


1 


1 


\p 2 - Pi + £ 3 ] 


2^/m\m 2 mi 

7 [~P2 + Pi + 6] , 

\Jm\m 2 

m 2 1 


mi \Jm\m 2 

2 £a 


\_P2 - Pi + £ 3 ] 


y/mim 2 


(74) 

(75) 

(76) 

(77) 
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The spectral densities of the form factor combinations in the l.h.s. of eqs. (| E 0 ) read 


9~s = - /3, 

m i 

9i + fi = -4/3, 
ff 2 + /2 = 4 ^/ 3 , 

mi 

Comparison with the eqs. dzl © yield the relation 

PiM = p 2 (u), 

For the form factor £ 3 we come to the representation 


1 

/ [‘ djj 1 

dz 2 Mz2)V^T2^) / -yMz!)- 


m? — 


— (m 3 + zi)(m 3 + z 2 ) + ^ 

1 UJ A — 1 


(78) 

(79) 

(80) 

(81) 

(82) 


with z\ given by the expression (|33|). In ( p2| ) we have neglected the 0(1 /viq) difference between oj and u 
On the other hand, the equation ([77|) yields the representation for the form factor £3 in a different form 


1 

6 M = J dz 2 fa(z 2 ) y/z 2 (z 2 + 2 m 3 ) J ^</>o(^i) 


z 2 + 2 m 3 

u + 1 


z 2 + m 3 + 


(zi - z 2 ) ui 

UJ 1 


m 3 

“ 2 “ 


to 3 


2 m 3 + 2:1 + 22 
1 + u> 


(83) 


One can check that for the soft wave functions providing convergency of the integrals and nonsingular at 2 = 0 
the representations (^2|) and (|83|) are equivalent. At zero recoil one finds 


6 ( 1 ) = J dz 2 \/z 2 (z 2 + 2 m 3 ) 3 / 2 <fe)f = 


(84) 


where ( 2 ) is an average kinetic energy of the light constituent quark inside an (infinitely) heavy meson in its 
rest frame. It is worth noting that the function £3 is positive for all uj. 

The universal form factor pi = p 2 can be found from the expansion of s(u>) with a NLO accuracy. In this 
case the NLO terms in the 1 /toq expansions of the integration measure, the wave function, and the spectral 
density s contribute. Namely, we can write for the form factor h s the expression 


h s (uj) = J 


dp 0 


dp[ 1} 


dp[ 2) 

m2 


S, S 2 

s[uj) H-1- 

m 1 m 2 


M z i) ( 1 + £^Xi(*i) ) M z 2) 


1 + ^-Xx (z 2 ) 


4 m 2 ' 


= + 


— + —) e (1) + = + K (1) M + A(w - 1)] 

TO 1 TO 2 / \TOi TO 2 y 


(85) 


and hence pi(w) = £^(u;) + £'(w)A(u; — 1 ) — £ 3 (w) — We do not present explicit expression for pi. 

However, let us consider pi at zero recoil. The analysis of the anomalous contribution at uj = 1 gives 


h s (u) 





( 86 ) 


Using the relations h s (u> = 1) = h s (u> = 1) + 0(1/toq) and £ 3 (1) = ( 2 )/3 we obtain pi(l) = 0 just as required 
by the HQET. 


F. The NLO analysis of the form factor / 

First, let us demonstrate that the dispersion representation of the form factor / requires subtraction. To 
this end consider the anomalous contribution at u> = 1. For the form factor fu constructed from the spectral 
density fn through the unsubtracted double dispersion representation we find [] 


1 Hereafter we denote as /a the form factor constructed from the spectral density a through unsubtracted double 
dispersion representation. 
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foil) = 2y/M 1 M 2 


r, 1 

( 1 

1 \ 

i + - 


H-) 

2 

(mi 

m 2 J 


H-I (( 2 ) + m 3 — A) + ... 

TO 2 J 


(87) 


This value contradicts the Luke theorem which requires the 1 /rriQ corrections to vanish at zero recoil. Let us 
demonstrate that the form factor / constructed from the spectral density with the included subtraction term 

/ = Id + (Ml - Si + M| - s 2 )g = Jd + (2pip 2 - 2pip 2 )g, (88) 

satisfies the NLO relation ©• Here 2pip 2 = si + s 2 — S 3 . 

We have noticed above that to LO and NLO the expansion of the form factor / and the combination 2pip 2 ■ 
g — TOi ■ g 2 coincide (|2l]). Hence checking the NLO relations for the form factor / is equivalent to checking with 
the NLO accuracy the relation 

fo — f 2 p 1 p 2 g ~ mi ■ g 2 . (89) 

The spectral density of the r.h.s. of eq. (89) can be written as 

2p\p 2 ' g mi • g 2 — 2m\m 2 (u) + l)s + [ 2 mi(z 2 + m 3 ) + 2m 2 (zi + m 3 )]s - 8 (mi + m 2 )/3 — 4 m 2 {u> — 1 )/3. (90) 

in NLO we need the expansion of the spectral density s in LO and NLO which 


For checking the expression 
has the structure 
1 


-s = ni\cx 2 + m 2 a 1 + m 3 (l — ai — a 2 ) = m 3 + 


zi + z 2 + 2 ro 3 


id +1 


s L+ s 2 _ 

mi m 2 


with 


Si = \z 2 {z 2 + 2 m 3 ) - 2 iTO 3 - (z 2 + 2 to 3 )—!:-^ 

2 id — 1 


(zi + z 2 + 2 m 3 ) 2 


id + 1 

and S 2 is obtained from S 1 by replacing z\ and z 2 . The spectral density fn reads 

Zi + z 2 + 2 ?u 3 


(91) 


(92) 


f D = 4toito 2 (w + 1) m 3 + 


id - 


1 


+ 2 m 2 z 2 (z 2 + 2 m 3 ) + 2rmzi{zi + 2 m 3 ) — 4m 2 (w — 1 )/3. (93) 


Notice that for checking the NLO relation (|89j) between the form factors we do not need explicit expression for 
the integration measure in the NLO: in the LO the spectral densities are equal and hence the NLO contributions 
from the integration measure into both sides of the eq ©) are equal too. Finally, the eq. (|89|) is satisfied if the 
following relation is valid 


1 

f c /?j 

f dz 2 \/z 2 (z 2 + 2 m 3 )(/>o(z 2 ) / -y(f>o(zi) 

-1 


-z 2 (z 2 + 2 m3) + 2 (zi + m 3 ) ( m 3 + 


z\ + z 2 + 2 m 3 
id + 1 


— 4/3 + 2(u) + 1) aS*! 


= 0 , 


(94) 


with Z\ given by (J63|). One can check this relation to be true for any function <p 0 (z) regular at z = 0. Hence, 
the form factor / calculated with the subtracted double spectral density ( |S 8 |) satisfies the HQET relations in 
LO and NLO at all id. 

Strictly speaking, the NLO analysis does not allow us to uniquely specify the subtraction term: namely, any 
spectral density of the form 

/ = Id + (Mj 2 — Si + Ml — s 2 )pf (95) 

has a proper NLO behavior in accordance with © provided the spectral density pf behaves in the LO as 

„ / zi + z 2 + 2 m 3 ' 

« = 2 ( m3 0+1 

As we demonstrate in the next section, the analysis of the heavy-to-light LO relations requires the identification 

PS =9- .... n 

Let us notice that the form factor obtained within the light-cone approach 0 can be represented as dispersion 
representation with the spectral density [lq] 


he = ^=f D 

VS2 


Mo 


si - s 2 — s 3 Ml - Ml - s 3 \ ai + 


2a/S2 


2 M 2 


a 2 


At zero recoil one finds 

/lc( 1) = 2y/MiM 2 
that contradicts the Luke theorem. 


1 - - 

f± 

— ) 

2 

[mi 

m 2 ) 


+ m 3 - A) + ... 


(96) 


(97) 
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G. The NLO analysis of the form factors ai,2 and 


The unsubtracted spectral densities of the form factors ai and a' 2 = a 2 — 2s in the LO read 


di d = 4 [(22 + 2m 3 )ai - m 2 an(w - 1)], 


(98) 


® 2 D ~ 4 


TO 2 , - . 

— (22 + 2 m 3 )ai 
m 1 


m 2 ai 2 {io — 1) 


(99) 


The quantity 02 is more convenient than a 2 for calculations as the LO term of the heavy quark expansion of a' 2 
is zero. 

The unsubtracted spectral representation for ai(u>)m 2 — o 2 (w)mi in combination with eq (|23|) gives 


X 2 M = 


1 

1 /* fd 

- / dz 2 Mz2)VM^T^T) / -y</>o(2l) [ 81^2 ^ d 2 Wi] 

-1 


( 100 ) 


Substituting the representation (|||) for £ 3 we find 


X 2 M = 0. (101) 

Let us now consider the linear combination ai(uj)m 2 + a' 2 (u>)m\. As a first step, show that the unsubtracted 
dispersion representation is not compatible with HQET. To this end calculate the unsubtracted form factors 
air, and a' 2D at zero recoil: 


ai£>(l) = 

1 

1 

2 (a 

x + 7X13 

sjmim 2 

m 2 

|_3 

)+ 2 

a 2£>(l) = 

1 

1 

1 (- 

x + m3 

i/mim 2 

mi 

[3' 

2 


and hence 


aiD(l)m 2 + a 2D (l)mi = 


1 


v /mim 2 


\{z) + m 3 \ 


in contradiction with the HQET result eq. (p 


A 


y/rmm2 


( 102 ) 

(103) 

(104) 

(105) 


This fact suggests a necessity of subtraction in the quantity aim 2 + a' 2 m\. Let us write the spectral density 
with subtraction in the form 


aim 2 + a 2 m\ = aiDrn 2 + a 2D mi + — 


k f Mi — Si M 2 — s 2 \ p, 


with 


p a ~ 2 m 3 + 


w + 1 \ 2^/sT 


Z\ + z 2 + 2 to 3 


+ 


2 i/S 2 


(106) 


to + 1 

Then corresponding representation for the form factor reads 


_ _L 

ai(u;)m 2 + a' 2 (u)mi = v ^ im2 / dz 2 (/)o{z 2 )y/z 2 (z 2 + 2 ro 3 ) / 0 ( 21 ) 


1 ™.3 


-1+Z2+2TO3 \ 

u+1 I 


x a lfl m2 + a 2 D m 1 - -^{zi +22 + 2m 3 ) ( 

According to the HQET eq (^4|) this relation should be equal to 

1 r 2A u-1 


I kA g(A) 

' 2y/mirrL2 tj+l ' 


y/mim 2 \_u> + 1 


w + 1 


(107) 

(108) 

(109) 
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The term proportional £(w) yields k = 4. One can check that this value also makes other parts of both 
expressions equal. Hence we arrive at the subtracted spectral density 


diTO2 + a 2 mi = di£)TO2 + «2D m i + 


1 / Mi — si M| — S2 


Pa- 


( 110 ) 


tO + 1 V v^T \/S 2 

The resulting spectral densities of the form factors ai and 02 with the built-in subtraction terms take the form 

1 

01 = am + - - 

(to + 1 

1 

a 2 = o,2d + 


(to + l)mi V y/Sl 


Pa 

2 ’ 

(111) 

Pa 

2 ' 

(112) 


H. The NLO analysis of go. 


The unsubtracted spectral density <?oz> has the form 

4 


9od — 


y/mim,2 


m 2 ai(m 3 + m 2 ai + mia 2 ) - 


P 


to + 1 


At zero recoil one finds 


Sod( 1 ) = 


1 


(nnm 2 ) 3 / 2 [ 2 


(z) +m 3 (z) 

6 


On the other hand, taking into account our earlier finding \ 2 = 0, the HQET result (20) reads 

1 


5od(1) = 


(TO 1 TO 2) 3 / 2 


A (£) 
2 6 


(113) 


(114) 


(115) 


Our experience obtained in considering aim 2 + a' 2 m\ hints that the subtraction procedure adds a term pro¬ 
portional at zero recoil to A — (z) — m 3l and hence we expect subtraction to work properly also in the case of 
9o- 

As a matter of fact, the spectral density 


go = goD + z— 


Ml - si | M 2 __s 2 \ Pgo_ 

J 2 


(to + l)mim 2 V \fs 1 
satisfies the HQ expansion (^) provided the function p go behaves in the LO as 


(116) 


Pgo - 2 [m 3 + 


zi + z 2 + 2 m 3 
to +1 


Concluding this section let us summarize our main results: we have calculated the universal form factors and 
demonstrated the dispersion representations with relevant subtractions in the case of heavy-to-heavy transitions 
to reproduce the structure of the heavy-quark expansion in QCD in the leading and next-to-leading orders. 
However, we have not been able to fix uniquely these subtractions. As we shall see in the next section, the 
heavy-to-light transitions provide further restrictions on the form of the subtraction terms. 


V. HEAVY-TO-LIGHT MESON TRANSITIONS 

In this section we discuss meson decays induced by the heavy-to-light quark transitions in which case Mi = 
m 2 + 0(1) is large, while M 2 ~ mi is kept finite. As found by Isgur and Wise |j^], in the region to — 1 = 0(1) 
the form factors of the tensor current can be expressed through the form factors of the vector and axial-vector 
currents in the leading l/m 2 order as follows 
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(117) 

(118) 
(119) 


s(q 2 ) 
92(q 2 ) 
9 o(q 2 ) 

9 i(<l 2 ) 


2 mJ 1 ^ ^ 

-2 M l 9 (q 2 ), 

^[2s(g 2 ) + a 2 (q 2 )], 

^-[-/O? 2 ) + %PiP2 -g(q 2 ) 


( 120 ) 


We address the two issue s: (i) pe rform the leading order l/m^ expansion 
fulfillment of the relations ( 117 12C| ) and (ii) discuss the scaling behavior of 
of different heavy mesons into a fixed final light state. 


of the form factors and show the 
the form factors of the transition 


A. The LO I/uiq expansion of the form factors in the quark model 


In the case of heavy-light transitions one observes the appearance of the two scales: the light scale M 2 — 
mi — m 3 — A, and the heavy scale M 2 — mi, and we may expand the form factors in inverse powers of the 
small parameter A/m 2 . The kinematical coefficients in the leading A/to 2 order in the region Co — 1 = 0(1) 
simplify to 


A(si ,s 2 ,q 2 ) 
ai 


a 2 


Oi 12 


an 


4to2 [(21 +m 3 +wtoi) 2 - s 2 ] , P = 0{ 1), 

1 (zi + m 3 + (Dtoi)(s 2 - mf + m §)/2 — s 2 {z\ + m 3 ) 


m 2 (z\ + m 3 + Mi ) 2 — s 2 

(zi +m 3 + u>mi)(zi + m 3 ) — (s 2 — m\ + m §)/2 


= 0 


OL\OL2 — 


(z 1 + m 3 + Mi ) 2 — s 2 
(3(zi + m 3 + Comi) 

2 m 2 [{zi + m 3 + Mi) 2 - s 2 ] 
1 /3s 2 


= 0 ( 1 ), 


= 0 


1 

m 2 


m 2 [(~i + rn 3 + Mi ) 2 - S 2 ] 


= O 



( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 


For the double spectral densities /» these expressions yield the following LO relations 

1 


2771 ' 


■h, 


g 2 = 2 m 2 g, 

9od = —(2 g + o, 2 d), 
m 2 

9 i = —[-Id + 2p!p 2 ■ g\, 


(126) 

(127) 

(128) 

(129) 


with 2pip 2 = Si + S 2 — S 3 . Notice that these relations hold also for Co = 0(m 2 ). 

First two equations directly give the LO equality of the corresponding form factors. The relation (128) 
between the unsubtracted spectral densities is more interesting and requires the LO identity of the subtraction 
terms 


Pa — Pgo ■ 


(130) 


The choice p a = p go = g is acceptable al thou gh there are no firm backgrounds to justify this very choice. 

The most informative is the relation ( |l29| ): this relation not only suggests a necessity of subtraction in the 
form factor / but also determines the subtraction term. When considering the heavy-to-heavy transitions we 
have found that the subtracted spectral density of the form (| 8 g|) 


/ — Id + (Mi — Si + M2 — 82) pf 


matches the HQET expansion if pf 


g in LO. The eq. (|129| ) prescribes pj = g. 


(131) 
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B. The scaling of the form factors with ttiq 


Let us point out that in the case of the heavy-to-light transition there is no parametrical suppression of 
the anomalous contribution to the form factors compared with the normal one as it was in heavy-to-heavy 
transitions and both parts contribute on equal footings. We shall present all the results for the normal part but 
the same holds also for the anomalous part. 

In the region Co = 0(1) one finds for the normal contribution 


/»M = 


y/rn2 , 


ds 2 tp{s 2 )\ l/2 (s 2 ,ml,ml) 
167r mi 


dgVu ;' 2 - 1 <j>o{z 1 )fi(z 1 ,Z 2 ,m 1 ,m 2 ,m 3 ,uj) 
i/(m 1 (<D + 1) + Z\ + z 2 + 2 m 3 )(nii(u) — 1) + z\ — z 2 ) 


(132) 


where 


s 2 = (mi + m 3 + z 2 ) 2 , 

Zl = Z 2 CC (l + -- 3 + ^ + m 3 (Co — 1 ) + r/V Co 2 — l 

V 7771 J 


\ 1 / 2 (s 2 ,ml,m%) 

2 mi 


0 (l/m 2 ). 


Using the relations (121) one finds that the spectral densities scale at large m 2 as 

fi = m^Piiu), mi, m 3 , z\,z 2 ). 

Namely, 

fi = 0(1), f 2 = 0 (?t 7 2 ), g = 0( 1), ai = 0 (l/m 2 ), d 2 = 0 ( 1 ), 
/ = 0 (m 2 ), 3i=0(l), 3 2 = 0(?ti 2 ), 3 0 = O(l/m 2 ), s = 0 ( 1 ). 

Hence, the form factors have the scaling behaviour of the form 

fi = m 2 i- 1 / 2 ri(w;mi,m 3 ;<^ 2 ,((>o). 

The variables u and Co are connected with each other as follows 

M 2 A 

to = - io - f 0(l/m 2 ), 

m i m i 

and hence the ratio 


(133) 


(134) 


(135) 


fM/M 


m+ 1/2 


(136) 


is universal for the transition of any heavy meson into a fixed light state. This behaviour reproduces the results 
of up to the logarithmic corrections which arise from the anomalous scaling of the quark currents in QCD. 

To summarize, matching the LO l/m 2 relations between the form factors in the quark model to the corre¬ 
sponding relations of j32| allowed us to determine the subtraction term in the form factor / and equated to 
each other subtraction terms in the form factors ay 2 and go- 
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VI. NUMERICAL ANALYSIS 


In this section we apply the derived results to the model-dependent estimates of the universal form factors. 

1. First, let us study the dependence of the IW function on the parameters of the quark model. We adopt 
the exponential parametrization of the radial wave functions in the form w(k 2 ) ~ exp(— k 2 /2(3 2 ). For obtaining 
the LO wave function </>o we must take into account the relationship between the parameters z and k 2 which is 
found from the following equation 


\fs = \Jk 2 + itiq + \Jk 2 + m| = ttiq + m 3 + z. 

This yields the relation 

k 2 = z(z + 27773 ) + 0(1/toq). 

Using the eq. (jii]) we come to the following form of the LO wave function 


4> 0 {z) ~ 


+ m 3 ( z(z + 27773 ) \ 

-“H- 20l~) 


Z + 2?773 

where (3 3 is the LO harmonic oscillator parameter 

P{m Q ) = /3 0 + 0{l/m Q ). 


(137) 

(138) 


(139) 


(140) 


For calculations we need to specify (3o- The exact value of (3q is not known but extrapolating the known values 
of Bp, [3 b , and (3d* we find /3o ~ 0.4 both in the WSB Q and the ISGW2 || models. In the model of ref. 
fl4[] from the analysis of the axial-vector decay constants of pseudoscalar mesons the following approximate 
dependence has been proposed 

(3{m Q ) ~ 2.5 mQ ^ 13 , (141) 

mQ + m 3 

which gives (3o ~ 2.5 m 3 . Table | presents the relevant numerical parameters. The results of calculating the IW 
function through the eq. m are shown in Fig. |j. Table || demonstrates the values £'(1) calculated with the 
eq. © and the parameters of the quadratic fit of the form 

CM =l-p 2 i (u>-l) + Si(u- l) 2 , (142) 


obtained by interpolating the results of calculations in the range 1 < to < 1.5. One can observe the value of 
—£'(1) to be considerably larger than the parameter p 2 obtained by the interpolation procedure. 


TABLE I. Parameters of the quadratic fit to the IW function and the NLO form factor £3 in various quark-model 
versions. 


Ref. 

m 3 

Po 


P 2 

5 

6 (f) 

Pe, 

<^3 

Set 1 [: 


0.35 

0.4 

1.47 

1.34 

0.78 

0.077 

1.49 

0.88 

Set 2 jj 

L 

0.33 

0.4 

1.43 

1.31 

0.76 

0.08 

1.4 

0.84 

Set 3 [|T 

i 

0.25 

0.63 

1.37 

1.12 

0.62 

0.17 

1.2 

0.65 



(0 


FIG. 2. The Isgur-Wise function calculated with various quark-model parameters: dashed line - set 1, solid line - set 
2, dotted line - set 3. 
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The IW function obtained from the dispersion quark model agrees with the results of the light-cone quark 
model applied to the elastic P —> P transition for (infinitely) heavy mesons Q £(1.5) = 0.62 and turns out to 
be a bit smaller than the SR result £(1.5) = 0.66 |Q. 

The NLO form factor £ 3 ( 0 ;) is found to be very sensitive to the light-quark mass. The parameters of the 
quadratic fit to £ 3 ( 0 ;) calculated through eq. ( |32| ) in the form 

hi(u>) = hi{ 1) [1 - pj (w - 1) + 8 i{w - l) 2 ] (143) 

for various sets of the quark-model parameters are shown in Table |[ One can observe an approximate relation 
£ 3 (w)/£(w) ~ 0.08 GeV for the light quark mass TO 3 ~ 0.35 GeV and £ 3 (w)/£(w) — 0.17 GeV for the light quark 
mass m 3 ~ 0.25 GeV in the whole region 1 < u < 1.5. This ratio is to be compared with the SR result with 
the 0 (a s ) corrections omitted [ pl| 

£ 3 (w)/£(w) = A/3 ~ 0.16 GeV, A ~ 0.5 GeV. (144) 

2. We are in a position to estimate the higher order corrections as we can calculate the form factors at finite 
masses and the leading-order contribution separately. For the ISGW2 parameter set (m c = 1.82, mb = 5.2, 
Pd = 0.45, Pb = 0.43, and Pd- = 0.38), the results on hf + and hf are s hown in Fig. ||, and Table [fi] presents 
the results of interpolating in the range 1 < u> < 1.5 with a quadratic fit (|143|). 


TABLE II. The form factors calculated with the parameters of ISGW2 model. 



hf+ 

h f 

£ 

Mi) 

0.93 

0.96 

1.0 

p 2 

0.87 

1.06 

1.25 

<5 

0.37 

0.53 

0.7 



(0 


FIG. 3. The velocity-dependent form factors calculated for the ISGW2 parameters: 
dashed line - h/(u), dotted line - hf + (u>). 


solid line - the IW function, 


The form factor hf can be seen to agree with the results of a combined fit to the ALEPH, ARGUS, CLEO, 
DELPHI, and OPAL data @ 


p) = 1.07 ±0.20. 


(145) 


One observes sizeable difference between the form factors hf and hf + and the Isgur-Wise function both in the 
value at zero recoil and the slope parameter. This indicates sizeable higher order corrections. At zero recoil the 
1 /toq corrections to hf vanish and writing J 31 J 


h f (l) = l + S 1/m 2 (146) 

we find Si / m 2 = —0.04 that agrees well with the value — 0.055±0.025 |u]]. At the nonzero recoil the higher-order 
corrections yield a more fiat w-dependence of the form factor hf compared with the IW function such that the 
slope of the IW function in our model turns out to be considerably larger compared with the slope of ' hf . An 
approximate relation p 2 IW ~ p 2 f — 0.2 is found. This is different from the SR result p 2 IW ~ pf lf + 0.2 Q. Let 
us point out that h'( 1) turns out to be considerably larger than the parameter p 2 of the quadratic fit obtained 
by the interpolation over the range 1 < u> < 1.5 just as in the case of the Isgur-Wise function. 
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VII. CONCLUSION 


We have performed a detailed analysis of the dispersion formulation of the quark model for meson decays 
based on representing the transition form factors as double spectral integrals in the invariant masses of the 
initial and final qq pairs through the soft wave functions of the initial and final mesons. The unsubtracted 
double spectral densities of the form factors at spacelike momentum transfers are calculated from the Feynman 
graphs, and subtractions in the spectral representations are determined by considering the two following cases: 
(i) by performing the heavy-quark expansion of the form factors in the case of heavy-to-heavy transitions 
and matching them to the HQET and (ii) by considering the relations between the form factors of the vector, 
axial-vector and tensor currents in the case of the heavy-to-light transition and matching them to the general 
relations of ref. |^|. This procedure gives the spectral representations with appropriate subtractions at space¬ 
like momentum transfers. The form factors at timelike momentum transfers are obtained by performing the 
analytical continuation in q 2 . The analytical continuation yields the appearance of the anomalous contribution 
to the form factors at q 2 >0. As a result, we have come to an explicit model of the form factors describing 
the meson decays induced by the transition between finite-mass quarks which develops the correct structure of 
the heavy-quark expansion. Let us emphasize that the obtained representations allow us to calculate the form 
factors directly at the timelike region. 

Notice that for deriving the heavy-quark expansion we impose only a strong localization of the meson mo¬ 
mentum distribution with a width of order of the confinement scale. No other constraints on the soft wave 
functions or numerical parameters of the model have emerged. 

Our main results are: 

1. We apply the method of ref. [^9| to meson amplitudes of the tensor current and find the expansions of the 
relevant form factors in LO and NLO in HQET. The 1 /mg corrections in the form factor h g+ are found to 
vanish at zero recoil just as in hf + and hf. 

2. The 1 /toq expansion of the soft wave function is constructed. The wave function normalization condition 
which is a consequence of the vector current conservation in the full theory yields an infinite chain of normal¬ 
ization conditions for the wave function components emerging in various 1 /toq orders. The LO normalization 
condition which is connected with the current conservation in the effective theory provides the correct normal¬ 
ization of the Isgur-Wise function at zero recoil. 

3. In the LO all the transition form factors are represented through the Isgur-Wise function in accordance with 
HQET. The Isgur-Wise function is calculated through the LO component of the heavy meson wave function. 
The subtraction terms in spectral representations of the form factors do not contribute in LO. The anomalous 
contribution which emerges when the analytical continuation to the timelike region is performed comes into the 
game only in the 1 /rriQ vicinity of the of the zero recoil point but otherwise is suppressed by at least the second 
inverse power of the heavy-quark mass. 

4. The NLO analysis of the P —> P transitions shows that the unsubtracted dispersion representations develop 
1 /tuq structure in accordance with HQET and allows calculating £3 through the LO wave function. 

5. In the case of the P —> V transition the unsubtracted representations for the form factors g , g\, < 72 , and 
777.102 — 771201 have the NLO behavior compatible with HQET and yield the relation p\ = P2 (or X3 = 0) and 
X ‘2 = 0. At the same time, the unsubtracted form factors /, go, and mia 2 + 771201 do not agree with HQET in 
NLO. The matching procedure allows restricting the form of the subtraction terms. 

6 . Analyzing the heavy-to-light transitions and requiring the fulfillment of the Isgur-Wise relations between 
the form factors of the tensor and vector and axial-vector currents [ fi 2 | further constrains the subtraction terms. 

7. We observe a discrepancy between the predictions of the various versions of the quark model on the NLO 
universal form factors : Namely, the analysys of the WSB model [Q performed in Q resulted in pi ^ p 2 7 ^ 0 
and P3 ^ 0, but P4 = 0; the quasipotential quark model || predicts all the NLO form factors to be nonzero; a 
consistent relativistic treatment of only qq intermediate states in Feynman graphs performed in this paper gives 
Pi = P2 , P3=0, and p 4 ^ 0. The origin of this discrepancy between the quark models is not fully understood 
and should be considered in more detail. 

8. The numerical results of the dispersion quark model using several parameter sets (i.e. constituent quark 
masses and wave functions) demonstrate a moderate dependence of the Isgur-Wise functions on the light-quark 
mass and the shape of the wave function. Namely, the light-quark mass 777.3 = 0.25-^0.35 and the wave-function 
width {3q = 0.4 -7- 0.65 yield £(1.5) = 0.55 -j- 0.6 that is a bit smaller than the SR result £(1-5) = 0.66 |Q. 

The form factor hf is found to have the behavior in agreement with predictions of other models and experi¬ 
mental results. We observe sizeable difference between the form factor hf and the Isgur-Wise function both in 
the value at zero recoil and the slope parameter. 

The size of the higher-order 1/mQ-corrections to hf( 1), <5i/ m 2 = —0.04, obtained using the ISGW2 parameters 
agrees favorably with the sum rule estimates. 

The slope of the IW function in our model turns out to be considerably larger compared with the slope of 
hf. An approximate relation p 2 w ~ p 2 f — 0.2 is found. This is opposite to the SR result p 2 IW ~ pf lf + 0.2 
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|3l|| . We would like to notice that the slope parameter is very sensitive to the interpolation procedure: namely, 
the value h'( 1) turns out to be considerably larger than the result of the quadratic interpolation over the range 
1 < w < 1.5. 

The dispersion quark model for the transition form factors can be further improved by performing the heavy- 
quark expansion in higher orders and matching to the HQET order by order. It should be taken into account 
that the quark model effectively describes the whole 1 /rriQ series, but the short-distance corrections are not 
contained in the model. The inclusion of such short-distance corrections into consideration is necessary for the 
application of the model to the analysis of the experimental results. 
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